A Blume-Capel Ising Ferromagnet with Annealed Vacancies on a Hierarchical Lattice 
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A dilute Ising ferromagnet is considered in this study using renormalization group techniques with 
a hierarchical lattice. A series of phase diagrams have been produced that probe the effects of varying 
the temperature and concentration of nonmagnetic impurities. Each phase diagram corresponds to 
a different strength for the internal coupling coefficients on our lattice. Phases have been interpreted 
and critical exponents calculated for the higher order transitions. 
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I. INTRODUCTION 

The Blume-Capel model is a spin-1 Ising model with a 
Hamiltonian having bilinear (J^ and crystal- field (Ay) 
interactions. 



with = 0,±1 (1) 

The concentration of nonmagnetic impurities (s; — 0) , 
or annealed vacancies, in our system is directly related 
to the crystal field interaction (~ A/ J). While the tra- 
ditional Ising interaction (Jij) primarily effects magnetic 
ordering. Systems driven by fluctuations in both magne- 
tization and density are particularly well suited for study 
using this model. Each contribution to the Hamiltonian 
in Eq. 1 involves a summation over nearest-neighbor (ij) 
pairs of our sites on our lattice unit structure including 
the crystal-field interaction term. Spin-1 Ising models, 
with density as an additional degree of freedom, have 
been used to probe, investigate and further understand a 
range of very complex systems. Included amongst these 
is the superfluid transitions in He^ — He"^ mixtures 
structural glasses [1], binary fluids, materials with mo- 
bile defects, semiconductor alloys [3] , microemulsions [3| , 
frustrated percolation systems [5] and aerogels @. 

Many different types of competing interactions have 
been the focus of previous studies using the Blume- 
Emery-Grifhths model in conjunction with mean-field 
methods M. I sl. llOl and/ or renormalization-group tech- 
niques prl [l3.1l3l. Il4 fill. [TgI. [TtI . In many studies, it is 
found that competition between underlying microscopic 
interactions can drastically alter phase diagrams and as- 
sociated criticality in various Ising systems. 

Competing bilinear interactions [l^l in a spin- 1/2 Ising 
model, competing bilinear interactions in a BEG sys- 
tem 16], competin g bi quadratic interactions in a dilute 
Ising ferromagnet [l7[, and simultaneous competition 
between crystal-field and biquadratic interactions in a 
BEG ferromagnet [15] have all been investigated using 
renormalization-group techniques in concert with hier- 
achical lattices. 



Other studies have considered the effects of quenched 
random bonds [3] and quenched random fields [l^ 
upon the criticality and phase diagrams in BEG systems. 
Branco et al. [Ij, [Ij, [20] have considered the effects of 
random crystal fields using real-space renormalization- 
group methods and mean-field approximations for both 
Blume-Emery-Griffiths and Blume-Capel hamiltonians, 
respectively. 

The current study complements these earlier works as 
it considers a dilute Ising ferromagnet and the effect of 
varying the concentration of nonmagnetic impurities in 
the system. A series of phase diagrams have been pro- 
duced while varying the relative strength of internal cou- 
pling coefficients in the system. For the model considered 
(dilute Ising on a hierarchical lattice) and techniques em- 
ployed (renormalization group theory) the existing liter- 
ature lacked data, thus, the purpose of the results and 
discussion presented here. 



II. RENORMALIZATION GROUP AND 
HIERARCHICAL LATTICES AND 

In general, infinite hierarchical lattices are constructed 
from a basic unit, or generator, by repeatedly replacing 
each bond by the basic unit itself. A generic hierarchical 
lattice is constructed and illustrated in Figures 1. Fig- 
ure 2 illustrates the construction of the more complex 
hierarchical lattice [2l|, used for the present study. 
The increased complexity due to changes in the inter- 
nal connectivity of the basic generating unit structure. 
Since these specialized lattices yield exact renormaliza- 
tion group recursion relations, hierarchical lattices are 
very attractive to use as model systems. Thus, critical 
scaling exponents and phase diagrams can be calculated 
very precisely. Several previous studies have emplwed hi- 
erarchical lattices to effectively probe s pin glass [ll, , 
frustrated [isl . [TtI . random-bond [23], random-field 
[IB], directed-path [2(| and dynamic scahng [13] sys- 
tems. 

Renormalization reverses the construction process for 
the infinite hierarchical lattice as internal degrees of free- 
dom are eliminated by summing over all configurations of 
the innermost spin sites (represented by solid black dots 
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FIG. 1: An infinite hierarchical lattice generated from a basic 
unit (Berker and Ostlund [13]). 




Cr = 5^exp[-/3ff'] = ^i?K^^-,s;) (6) 
with Ri{s^^ Sj) — exp[J'5^5j 

S ■ ,S ■ 

+if's?4-A'(.2 + ,2) + G'] (7) 

where G' is a constant used to calculate the free energy. 

The derivation for the renormalzation-group trans- 
formation for each coupling coefficient is calculated 
by equating individual contributions, Ri{si,Sj) and 
Ri>{si, Sj), to the summation for the partition function 
at each length scale. These contributions to the par- 
tition function, Ri(si,Sj) and Rii{si, Sj), correspond to 
the same fixed configuration of end spins, Si,Sj, at the 
two different length scales, I and V . From the resulting 
relationships, we algebraically derive relations between 
interaction strengths at the two length scales, I and 
J' (J, K, A), J' (J, K, A), A' (J, A). The reader is directed 
to Section 5 for a derivation of these relations. 

Phase diagrams are mapped and the order of each tran- 
sition determined using these recursion relations in con- 
junction with the initial values of J and A, and the re- 
sulting flow and sink(s) of the renormalization-group tra- 
jectories. 



FIG. 2: Construction of the hierarchical lattice used in this 
study. (Reprinted from Journal of Magnetism and Magnetic 
Materials, 314, 69-74 (2007) D. P. Snowman, with permission 
from Elsevier) 



in Figures 2a and 2b, represented by Si,Sj in Equation 

5). _ 

Renormalization-group relations, relating the coupling 
coefficients at the two length scales, are developed by re- 
quiring the partition function to remain unchanged with 
each rescaling. The new effective coupling coefficients J', 
and A' are separated by a distance V which is h lattice 
constants in the original system, where b is the length 
rescaling factor of the renormalization-group transfor- 
mation. Under renormalization, biquadratic interactions 
arise that must also be considered. Thus, in general. 



J'^Rj{J,K,^) (8) 
K'^Rk{J,K,A) (9) 
A' = Ra{J,K,A) (10) 



With each phase there exists an associated phase 
sink, see Table I, at which the values of the interac- 
tions (J, K, A) have reached a fixed point denoted by 
{J* , K* , A*). At these fixed points the system is scale 
invariant and as a consequence renormalization does not 
effect the properties of the system as the length scale is 
increased by a factor of b. That is, the fixed points must 
satisfy the recursion relations such that 

J* ^ Rj{J*,K*,A*) (11) 



Ci,{J',K',A')^CiiJ,K,A) 



(2) 



with /' = bl 



(3) 



K* ^ Rk{J*,K*,A*) (12) 



Ci = J2 (^M-m = E ^'(sm s,) (4) 



'3' - 

<ij> 



A* ^ RAiJ* , K* , A*) (13) 



The hierarchical lattice used in the current study 
with i?;(si, Sj) = exp[-/3i7] (5) is generated from a unit structure with two types of 

components (see Fig. 2a and b), similar to references 
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[H, m, H [13, nil. One type of component, type A, is 
distinguished by a cross-Hnk feature of strength p thus 
allowing internal spins to interact via nearest neighbor 
interaction J, A. The limitin g ca se oi p = corresponds 
to the hierarchical equivalent [21[ of the Migdal-Kadanoff 
[2^ Is^ l decimation-bond moving scheme in two dimen- 
sions. 

In addition, a second type of component, type B as 
shown in Figure 2b, allows the end spins to interact via 
two different connecting paths consisting of toi and m2 
pairs of spins, respectively. The relative number of each 
type of component, used in our basic generating unit, is 
controlled via two parameters pA and pb- 

III. PHASE TRANSITION 
CHARACTERIZATION 

The free energy density (dimensionless free energy per 
bond), /, can be expressed as 



where F is the Helmhotz free energy and Nf, denotes 
the total number of bonds in the system. The free en- 
ergy density consists of a sum, over all iterations of the 
renormalization-group transformation, of the contribu- 
tions G'^"-' to the free energy density due to the degrees 
of freedom removed during each transformation. Each 



renormalization-group transformation reduces the length 
scale of the system by a factor of b and the number of 
spins by a factor of b'^. 

Numerically differentiating the free energy density al- 
lows us to calculate densities, magnetizations, and near- 
est neighbor correlations. For example, the magnetiza- 
tion, m = = ^ jjj- , can be calculated by numerically 
measuring the shift in the free energy density with a small 
perturbation in the magnetic field, where Ns is the num- 
ber of sites. Similarly, the density can be calculated by 
differentiating the free energy density with respect to the 
crystal field coefficient, P = • Nearest neighbor 

correlations of the bilinear, (siSj) = ■^Ij, exchange in- 
teractions are also valuable when interpreting the phases 
and characterizing transitions. 

The order of each transition is characterized using the 
four thermodynamic quantities discussed above with first 
order transitions being signaled with discontinuities in 
any one of these order parameters. Note, second order 
or continuous transitions exhibit no such discontinuities. 
In addition, critical scaling exponents can be calculated 
at critical transitions since we have exact recursion rela- 
tions. 

IV. RECURSION RELATIONS 

For each fixed end-spin configuration, we equate the 
contributions to the partition function from the two 
length scales, allowing us to write the following equal- 
ities for the type A structure shown in Figure 2a, for the 
general case with nonzero K . 



= exp[-4A] +2exp[-2J + 2i^- A(6+p)] 

+2 exp[2 J + 2K^ A(6 + p)] + cxp[J(-4 + p) + K{4 + p) ^ A(8 + 2p)] 

+2 exp[- Jp + KiA + p)- A(8 + 2p)] + exp[J(4 + p) + K{A + p) - A(8 + 2p)] 

= exp[J' + - 2A' + G] =Rv[l,l], (15) 



Ri[l,Q\ = cxp[-2A] + 2cxp[-J + ii:- A(4+p)] + 2cxp[J + A'- A(4 + p)] 

+ cxp[J(-2 +p)+ K{2 +p)- A(6 + 2p)] + 2 exp[- Jp + K{2 + p) - A(6 + 2p)] 

+ exp[J{2+p)+K{2+p)-A{6 + 2p)]=cxp[-A' + G] = Ri,[l,0], (16) 

i?,[l,-l] = exp[-4:A] + Aexp[2K - A{6 + p)] 

+2 exp[- Jp + K{A + p) - A(8 + 2p)] + 2 exp[Jp + K{A + p) - A(8 + 2p)] 

= exp[- J' + K' - 2A' + G] = Ri- [1, -1], (17) 

i?i[0,0] = l + 4cxp[-A(2+p)] + 2cxp[~Jp + is:p- A(4 + 2p)] 

+2 exp[Jp + Kp-A{'i + 2p)] = exp[G] = R^ [0, 0] , (18) 
I 



The desired renormalization-group transformations, scales for the type A unit structure, can be derived by 
relating the coupling coefficients between the two length 
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algebraically manipulating the relationships above (Eqs. 
15-18). 



Ja = t; log 



2 ^i?p(l,-l) 
1 Ri,{l,l)Ri,il,-l)Rf,{0,0) 



= log 



R.f,{hO) 
Ri' (0,0) 



Ri' (1,0) 



(19) 

(20) 

(21) 
(22) 



The recursion relations for the less complex (type B) 

unit structures have the same form as in Eqs. 19-22, 
but the expressions (Eqs. 15-18) for the corresponding 
Rl{si,Sj) differ. Combining the contributions from both 
types of unit structures (type A and type B as shown in 
Fig. 2), the renormalization relationships become 



J = paJa+PbJb 
k' = paK'^+ pbK'b 
A' =paA'a+PbA'b 



(23) 
(24) 
(25) 



The exact nature of the rcnormalization-group rela- 
tions above, allow us to calculate critical exponents by 
linearizing the recursion relations in the vicinity of the 
critical transition under investigation. That is, 



J - J* 



K -K* 



Tjj{J-J*)+Tjk{K-K*) 
+Tja(A - A*), (26) 

TKj{J-r)+TKK{K-K*) 
+Tka(A-A*), (27) 



A - A* = 



Taj{J -r)+TAK{K-K*) 
+Taa(A- A*), 



(28) 

where Tjj = Taj = etc. and are evaluated 
at the fixed point in question. A recursion matrix, with 
elements Txy and eigenvalues of the form 



A, = 6^' (29) 

can be used to represent the critical relations in Eqs. 
26-28. Here, h is the length rescaling factor (in our case 
6 = 2) and yi represents the corresponding critical ex- 
ponent for the Z*^ eigenvalue. Critical scaling exponents 
have also been calculated for nonzero perturbations of 
the odd sector components H and L. 



Phase 


Sink 


Characteristics 


Dense Paramagnetic 


J ^ 





Low concentration of 




A ^ 


— oo 


nonmagnetic impurities 


Ferromagnetic 


J -> 


+ 00 


Magnetically 




A ^ 


—oo 


ordered 


Dilute Paramagnetic 


J ^ 





High concentration of 




A ^ 


-hoo 


nonmagnetic impurities 



TABLE I: Phases and Corresponding Sinks 



V. RESULTS 

Our results below probe the effects upon ordering of 
varying the temperature (1/J) and vacancy concentra- 
tion ( A/ J) present in the system. A series of phase 
diagrams are produced, each corresponding to a different 
level of internal interactions, timed via the parameter p. 
Exhaustive analysis of renormalization-group trajectories 
and corresponding sinks yields three phases: a ferromag- 
netic, a dense paramagnetic and a dilute paramagnetic 
phase. Each phase shares a common sink or basin of 
attraction in parameter space, as detailed in Table I. 

The dense and dilute paramagnetic phases are distin- 
guished from one another via the flow of the crystal-fleld 
interaction term. A flow to -oo corresponding to a dense 
population of magnetic species, whereas a flow to -|-oo 
corresponds to a system to dilute for magnetic ordering 
to occur. For each paramagnetic phase, the bilinear (J) 
interaction flows to zero indicating no magnetic order- 
ing. In addition, a ferromagnetic phase arises with a 
renormalization group flow with the bilinear interaction 
flowing to +00 and the crystal-fleld interaction flowing to 
-oo. The ferromagnetic phase can be reached by decreas- 
ing the temperature, or, by decreasing the concentration 
of nonmagnetic impurities. 

The results presented below were generated while con- 
sidering various planes of parameter space with different 
degrees of internal connectivity. The parameter p al- 
lows us to vary this c;onnectivity as we investigate the 
effects upon the underlying phase diagrams of changes in 
temperature (1/J) and crystal-field interaction ( A/J). 
In the first plane we consider the plane with p=l and 
we find two paramagnetic phases (dense and dilute) and 
a magnetically ordered ferromagnetic phase. The fer- 
romagnetic phase occurs at intermediate and lower tem- 
peratures for values of the crystal-field interaction (A/J) 
that correspond to a lower concentration of nonmagnetic 
impurities on our lattice. 

The ferromagnetic phase disorders at high tempera- 
ture via a second order transitions to a dense paramag- 
netic state. In this state, the concentration of occupied 
sites is such that magnetic order is possible, however, the 
temperature is too great for long-range magnetic order 
to propagate. The dense paramagnetic phase is present 
only at high temperatures and is separated from its dilute 
paramagnetic counterpart via a line of first order phase 
boundary that terminates at a high temperature critical 
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FIG. 3: Phase diagram with connectivity values 
{p,mi,m2,pA,PB) = (1, 8, 9, 40, 1), showing different 
basins of attraction and associated phases with critical 
endpoint (E) and critical point (C). Solid lines represent 
second-order transitions, whereas dashed lines represent 
first-order transitions. 



FIG. 5: Phase diagram with connectivity values 
{p,mi,m2,PA,PB) = (4, 8, 9, 40, 1), showing different 
basins of attraction and associated phases with critical 
endpoint (E). Solid lines represent second-order transitions, 
whereas dashed lines represent first-order transitions. 



First order pliase 
boundar-y tenninates 
at a liigh temperature 
critical point C 




-3.0 -1.0 1.0 
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5.0 



FIG. 4: Phase diagram with connectivity values 
(p, mi, m2,PA,Pfl) = (2, 8, 9, 40, 1), showing different 
basins of attraction and associated phases with critical 
endpoint (E) and critical point (C). Solid lines represent 
second-order transitions, whereas dashed lines represent 
first-order transitions. 



point C. At temperatures above this critical point, it is 
possible to drive the system from one paramagnetic phase 
to the other sans transition as in the standard liquid-gas 
phase diagram. The line of criticality separating the fer- 
romagnetic and dense paramagnetic states terminates at 
a critical endpoint E upon the intersection with the line 
of first order transitions. 

The dilute paramagnetic state appears at those val- 
ues of the crystal-field interaction that correspond to a 
concentration of nonmagnetic impurities that exceeds a 
certain percolation threshold necessary in order for long- 
range magnetic order to propagate. 

An increase in the strength of the internal coupling 
coefficients to p=2 reveals an underlying phase diagram 
that is very similar to the p=l case. The location of 
the first order phase boundary remains unchanged, and 
nearly vertical at A/ J — 0. Thus, our system can be 
forced to transition between the two paramagnetic states 
at high temperature only through a change in the crystal- 
field interactions present. That is, we must change the 
concentration of occupied sites to traverse this first-order 
phase boundary. Also, the line of criticality separating 
the ferromagnetic and dense paramagnetic phases has 
shifted to higher temperatures. 

In our final two phase diagrams we consider the ef- 
fects of increasing the internal connectivity to p=4 (Fig- 
ure 5) and p=8 (Figure 6). Topologically the underlying 
phase diagrams remain unchanged with this increase in 
the strength of the internal coupling coefficients. In each 
phase diagram investigated, however, there is evidence 
(at low temperatures) of an ordering transition driven 
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FIG. 6: Phase diagram with connectivity values 
{p,mi,m2,pA,PB) = (8, 8, 9, 40, 1), showing different 
basins of attraction and associated phases with critical 
endpoint (E) and critical point (C). Solid lines represent 
second-order transitions, whereas dashed lines represent 
first-order transitions. 



by a decrease in temperature that would result in the 
dilute paramagnetic phase ordering to the ferromagnetic 
phase via a first order transition over a very small range 
of crystal-field interactions. The shift in the second-order 
phase boundary separating the dense paramagnetic and 
ferromagnetic phases continues to be pushed to higher 
temperatures with increasing internal connectivity. 

In each plane in parameter space considered in this 
study, lines of critical transitions have been probed. Lin- 
earization of the recursion relations while maintaining 
scaling fields associated with J, K, A, and, nonzero per- 
turbations in odd sector contributions H and L, results 
in a recursion matrix. Calculation of the eigenvalues al- 
low us to extract critical scaling exponents (as discussed 
in further detail in Seciton 4. This analysis has been 
conducted for the line of criticality separating the fer- 
romagnetic and dense paramagnetic states at high tem- 
peratures for each degree of internal connectivity (i.e. 
p = 1,2,4, and 8). 

Associated with scaling fields J, K, A for the case of 
p = 1, we find two relevant eigenvalues A2 = 13.48 and 
A4 = 2.00; corresponding to critical scaling exponents 
of 2/2 = 3.75 and j/4 — 1.00, respectively; and, an irrel- 
evant eigenvalue with Ag = —1.48. Associated with the 
odd sector scaling fields, H and L, we find two relevant 
eigenvalues Ai — 5.00 and A3 = 2.00; corresponding to 
critical scaling exponents of yi — 2.32 and 2/3 = 1.00. 

For internal connectivity p = 2, associated with scal- 
ing fields J,K,A, we find two relevant eigenvalues A2 — 
15.54 and A4 = 2.00; corresponding to critical scaling ex- 



ponents of j/2 = 3.96 and 1/4 = 1.00, respectively; and, an 
irrelevant eigenvalue with Ag — —1.54. Associated with 
H and L we find two relevant eigenvalues Ai — 6.00 and 
A3 — 2.00; corresponding to critical scaling exponents of 
yi = 2.58 and ys = 1.00. 

For the case of p = 4, associated with J, K, and A, we 
find two relevant eigenvalues A2 = 19.63 and A4 = 2.00; 
corresponding to critical scaling exponents of y2 = 4.29 
and ?;4 = 1.00, respectively; and, an irrelevant eigen- 
value with Ag — —1.63. Associated with H and L we 
find Ai = 8.00 and A3 — 2.00; corresponding to critical 
scaling exponents of yi = 3.00 and 7/3 = 1.00. 

In the last plane in parameter space, corresponding to 
p — 8, associated with scaling fields J, K, A, we find two 
relevant eigenvalues A2 = 27.73 and A4 = 2.00; corre- 
sponding to critical scaling exponents of y2 = 4.79 and 
j/4 = 1.00, respectively; and, an irrelevant eigenvalue with 
Ag ~ —1.73. Associated with the odd sector fields, H 
and L, we find two relevant eigenvalues Ai = 12.00 and 
A3 — 2.00; corresponding to critical scaling exponents of 
yi = 3.58 and ^3 = 1.00. 



VI. SUMMARY 

In summary, this study considers a Blume-Capel Ising 
ferromagnet using renormalization group methods and a 
hierarchical lattice. Exact recursion relations are devel- 
oped, phase diagrams calculated and critical exponents 
extracted as the concentration of annealed vacancies and 
temperature are varied. 

The results presented here detail our investigation 
into the effects of varying the strength of the inter- 
nal coupling coefficients via the connectivity parameter 
p. Three unique regions were found in each plane con- 
sidered: dense paramagnetic, dilute paramagnetic and 
ferromagnetic. The high temperature dense paramag- 
netic/ferromagnetic boundary was found to be second- 
order. This critical line terminating at the critical end- 
point E. The qualitative topology of each underlying 
phase diagram remained unchanged with increasing in- 
ternal connectivity p. The only subtle change being the 
gradual increase, with increasing p, of the critical line 
separating the high temperature dense paramagnetic and 
ferromagnetic phases. 

The primary goal of this work has been, and remains, 
to develop a complete and better understanding of the 
effects of varying temperature and density upon order- 
ing and criticality in ferromagnets with a range of in- 
ternal coupling coefficient strengths. Previous studies 
by this author [l^, [l^, [l^ have considered the role of 
various types of competing interactions using the same 
renormalization-group approach on hierarchical lattices 
with similar connectivity parameters. Thus, the present 
study - complete with uniform ferromagnetic and crystal- 
field interactions - can also be used as a standard for 
comparison when probing the effects due to competing 
bilinear and/or crystal-field interactions. 
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